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I. GENERAL FRAMEWORK 


Neutrinos are elementary particles with masses and mixings. The presence of mixings and non-degenerate masses 
produces peculiar phenomena when these weakly interacting particles propagate in vacuum, in astrophysical or cos¬ 
mological environments. In vacuum the mass eigenstates, solution of the Dirac equation, interfere and generate 
oscillations of the electron, muon and tau flavor states, a quantum-mechanical phenomenon analogous to Kq — Kq 
oscillations [I]. In astrophysical and cosmological environments, flavor conversion phenomena take place that are fun¬ 
damentally different from vacuum oscillations. They originate from neutrino interactions with the particles composing 
the medium, including neutrino self-interactions. In exploding stars, flavor conversion is produced by steep variations 
of the matter density due to front and reverse shocks, or by matter density fluctuations associated with turbulence. 
Interesting theoretically, flavor changes in matter impact direct observations of neutrinos from astrophysical sources. 
Moreover such modifications are connected to key open issues. In particular, they are relevant for assessing the site(s) 
where the heavy elements are made, for unravelling the explosion mechanism of core-collapse massive stars or for de¬ 
termining the impact of sterile neutrinos in the universe. These facts explain why mechanisms and precise conditions 
for flavor modification are under close scrutiny. 

To this goal, quantum kinetic frameworks are needed that can properly account for neutrino flavor conversion in 
media. In fact, the presence of neutrino mixings often requires special care in using known theoretical methods to 
derive neutrino evolution equations. Approaches differ for their potential to give evolution equations at an increased 
degree of complexity and for their easiness of being applied in situations of observational interest. Available equations 
can be classified in two categories depending on the inclusion of terms proportional to the first power of the Fermi 
coupling constant Gp, or to the second (G|’). The former account, in particular, for coherent forward scattering, 
while the latter for incoherent collisions. The two types are: i) mean-field and extended mean-held equations ; ii) 
Boltzmann equations. Note, however, that the Boltzmann equations used for neutrinos often do not include terms 
proportional to Gp', in particular, in the transport equations used in core-collapse supernova simulations. The use of 
i) or ii) depends on how dilute is the system in which neutrinos are travelling. Applications of interest exist for the 
Sun, for core-collapse supernovae - 0-Ne-Mg, iron-core supernovae, collapsars - and for the early Universe. Studies 
can also be found in neutron stars, in neutron star-neutron star or black hole-black hole binaries. Since stellar layers 
are dilute, the mean-held approximation has always been employed in the study of neutrino havor conversion. 

In the Sun, the density is low and the mean-held approximation is sufficient to implement the neutrino interactions 
with protons, electrons and neutrons. The result is an effective mean-held Hamiltonian that depends, in particular, 
on the electron number densities. The inclusion of such a term in the evolution equations gives rise, under appropriate 
conditions, to the Mikheev-Smirnov-Wolfenstein effect HEl. This is the efficient conversion of solar electron neutrinos 
into muon and tau neutrinos if neutrinos evolve adiabatically through a resonance. The resonance and the adiabaticity 
conditions are fulhlled depending on the neutrino energy, the values of one of the mixing angles (historically called 
6 * 12 ), the values and sign of one of the mass squared differences and the background density and density gradient. It 
is experimentally established that the MSW effect produces a supplementary suppression of the high energy ®B solar 
neutrino hux, compared to the averaged vacuum oscillations that reduce solar neutrinos with energies below 2 MeV 
(see, e.g., [HIS]). Neutrino experiments have forged the MSW effect a reference phenomenon for studies of neutrino 
flavor conversion in media. 

0-Ne-Mg and iron-core supernovae shine almost all their gravitational energy as neutrinos of all flavors when 
they collapse. In a ten seconds’ burst typically 10®® neutrinos are emitted with an average energy of about 10 MeV. 
Supernova neutrinos are produced in the dense region within the neutrino-sphere, which is opaque to neutrinos because 
of its high density. Since the timescale between collisions is short, out-of-equilibrium effects need to be considered. In 
this regime, neutrino propagation is determined through the solution of (relativistic) Boltzmann neutrino transport 
equation (see, e.g.. Ref. 0 ) in two and, recently, three dimensions. Such supernova simulations are computationally 
demanding and have not included neutrino mixings so far. This was physically motivated by the existence of a 
separation of length scales, namely the mean free path versus the flavour change length scale. 
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Neutrino mixings effects take place in the dilute stellar outer layers under the form of MS^^Q resonances [7]. 
Their location is fixed by the measured large neutrino mass squared differences, i.e. = 7.6 10“^ eV^ and 

Amf 2 = 2.4 10“^ eV^, the supernova density and neutrino energies in the tens of MeV energy range. However 
neutrino self-interaction turns out to be relevant in these environment^ because of the large neutrino densities. 
They have drastically changed the picture: collective flavor stable modes and flavor instabilities have been shown to 
emerge outside the neutrino-sphere, where neutrinos start free streaming (see, e.g.. Refs. [HIS]). So it appears at 
the moment that flavor modification takes place outside the gain region behind the shock, even in presence of the 
neutrino-neutrino interaction (see, e.g., cni) and does not impact the supernova dynamics. However, we cannot yet 
exclude that approximations made in available studies do not hinder instabilities in the gain region. Were this the 
case, the argument on the separation of scales would need to be seriously reconsidered. 

Indeed the transition region between the high and low density regimes is being actively discussed. So far the 
neutrino-sphere has been considered a sharp boundary between the Boltzmann-treated region and the mean-field 
description. In the former, supernova neutrinos acquire thermal distributions due to collisions. The corresponding 
fluxes are taken as initial conditions in flavor studies and evolved afterwards using coupled non-linear mean-field 
equations. Only about 10% of the neutrinos still undergo some collisions outside the neutrino-sphere. Although their 
fraction is small they might play a significant role in the competition between collisions and flavor evolution. 

The possible emergence of flavor instabilities in the region where collisions still matter require the solution of 
transport equations including neutrino mixings as well as the mean-field contributions (terms proportional to Gp). 
As we will see the quantum Boltzmann equations for neutrinos have been derived in Refs. [IIHIl]- Clearly it is 
desirable to solve such equations in some realistic cases to put current belief on a solid ground, or perhaps to discard 
it. In the meanwhile, preliminary studie^are of value. For example, a schematic study to investigate the possible role 
of collisions has consisted in implementing a small fraction of backward neutrino flux in the mean-field approximation. 
The outcome is that signihcant changes of the neutrino flavor can occur, compared to the case with only forward 
neutrino fluxes m- 

Besides collisions, other two-body correlations might be important in the transition region. In fact two novel kinds 
of two-point correlation functions naturally arise in the most general mean-field equations. These are associated with: 
i) non-zero neutrino masses min] or neutrino magnetic moments [18) : ii) neutrino-antineutrino pairing correlations 
[mis]. Note that mean-field contributions originating from the neutrino mass have first been pointed out in Ref. 
m where terms proportional to the ratio {m/E^)'^ are explicitly given (for the pairing mean-field), with m and 
being the neutrino mass and energy respectively. Ref. m has shown the presence of contributions of the order of 
(m/E^) and called spin coherence the neutrino-antineutrino mixing introduced by non-zero masses, for Majorana 
neutrinos. In particular the authors have shown that anisotropy of the medium is necessary to have spin coherence. 
Similar correlations produce neutrino-antineutrino mixing, in presence of the neutrino magnetic moment coupling to 
stellar magnetic fields [18] . As a result of neutrino-antineutrino mixing, flavor conversion between the neutrino and 
antineutrino sectors can take place. In Ref. m neutrino-antineutrino pairing correlation^ are first introduced that 
are formally analogous to the pairing correlations among electrons in the Bardeen-Cooper-Schrieffer (BCS) theory for 
supraconductivity, or to neutron-proton pairing in atomic nuclei. Although clearly different in nature from i), these 
correlations also introduce neutrino-antineutrino mixing and require anisotropy of the medium to be non-zero[17j. In 
Ref. [17j the most general evolution equations are derived for inhomogeneous and homogeneous media for both Dirac 
and Majorana neutrinos. They include both contributions from neutrino-antineutrino pairing correlations and from 
the neutrino mass (referred to as helicity coherence). 

The study of neutrino self-interaction effects with mean-field equations reveals a rich phenomenology of neutrino 
flavor conversion phenomena in dense media. The intrinsic many-body nature of the neutrino evolution in presence of 
neutrino-neutrino interactions was already emphasized in Ref. EH, where it was first pointed out that such interactions 


^ A third one named /ir resonance does not influence the neutrino fluxes. 

^ Note that their role was first pointed out in the context of the early universe. 

^ Note that the neutrino-sphere are energy and flavor dependent. 

In Ref. [15| the competition between the equilibration time due to collisions and the one due to flavor modification is studied in a 
schematic model of spins. Note however that the Boltzmann equations used are actually reduced to mean-field. Moreover, as the authors 
state it is unclear how much the model employed really tells us about a realistic supernova case. 

^ These correlations are usually discarded with the argument that they are expected to oscillate fast around zero |13l I2Q| . As we will 
discuss, this argument does not necessarily hold. 
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introduce a non-linear refractive index. Various collective flavor conversion effects have been identified. While most 
studies are realised in a core-collapse supernova ”set-up”, some of these findings are relevant for coUapsars or neutron 
star-neutron star (or black hole-black hole) binaries [22] . Note that interesting conversion effects can also take place in 
neutron stars (see, e.g.. Ref. [23j) The first studies are based on the simplihed bulb model jSj. In this model, the initial 
conditions at the neutrino-sphere correspond to neutrino emission with spherical symmetry and azimuthal symmetry 
along any radial direction. In this simplihed scheme and taking the ’’single-angle” approximatiorj^ a synchronization 
among flavor isospins is found, followed by a bipolar instability that can be seen as a havor or gyroscopic pendulum 
in havor space. At the end of bipolar oscillations and depending on the neutrino energy, neutrinos undergo either full 
or no conversion with a spectral swapping of their huxes above (or below) a critical energy (the spectral split). Such 
a phenomenon is understood either as a MSW phenomenon in a co-moving frame, or as a ’’magnetic resonance-like” 
phenomenon [24] . Moreover, different kinds of instabilities due to the neutrino self-interaction in supernovae are being 
uncovered using linearization. These break, e.g., some of the symmetries of the bulb model such as the azimuthal one 
[2S], showing the need for more realistic models (see also Ref. |26|b 

Clearly investigations of the equations of motion and of their solutions in simplihed models uncover general prop¬ 
erties of neutrino havor evolution in media. However simulations might be necessary under the geometry, convection 
and turbulence in realistic conditions to determine the impact on the r-process, or if instabilities occur behind the 
shock, or to clarify if typical timescales between collisions and havor evolution can compete. This is necessary for 
future observations of neutrinos from an (extra)galactic supernova or from explosions integrated over cosmological 
redshift (the dihuse supernova neutrino background), which might be possible with upcoming large scale neutrino 
detectors. 

In the early universe context, the situation is different. First, the collision rate among the relativistic particles is 
in competition with the expansion rate of the universe in the radiation dominated era, until freeze-out. Second, the 
hypothesis of homogeneity and isotropy of the plasma are well satished. Therefore, quantum Boltzmann equations for 
neutrinos are usually solved, under these assumptions, in baryogenesis and leptogenesis applications (see for example 
Refs. [23 [21]) or at the epoch of primordial light elements abundances (see Refs. [23110] for a review and references 
therein). The competition between collisions and flavor mixings is responsible in particular for an equilibration of the 
neutrino degeneracy parameter (or neutrino asymmetry) among the different neutrino flavors. In particular this has 
the interesting consequence that the tightest constrain on electron neutrinos applies to other flavors. Note however 
that, according to current observations, the neutrino degeneracy parameter is very small. Another application of 
interest is to quantify the extra radiation in the Universe in particular as sterile neutrinos. 

The goal of the present manuscript is to discuss theoretical approaches and the mean-held, extended mean-held 
and Boltzmann equations used in the description of neutrino propagation in media. To this aim we will highlight 
works that are representative of the methods employed so far. Moreover, we focus on derivations obtained for 
astrophysical applications where, in particular, the conditions of homogeneity and isotropy can be broken. In Section 
I, we summarise the density matrix approach and we present a standard method to obtain mean-held and extended 
mean-held equations. Then we focus on the Born-Bogoliubov-Green-Kirkwood-Yvon (BBGKY) framework. In Section 
B, we discuss a derivation based on the state-path integral approach. Section IB is focussed on the Keldysh or Glosed- 
Time-Path (”in-in”) formalism. Finally Section IV is our conclusion. 


II. THEORETICAL FRAMEWORKS 

In the last two decades various theoretical approaches have been employed to derive neutrino evolution equations 
for astrophysical and cosmological applications [I2HI11 [13 HHIEB ElMl [3TH35] . Several works use a formalism based 
on the density matrix. Its application to neutrinos was first introduced in Ref. |29) and widely used afterwards (see, 


® In the ’’single-angle” approximation one assumes that neutrinos are emitted at the neutrino-sphere with the same angle and have similar 
flavor histories. In the ’’multi-angle” approximation neutrinos emitted at the neutrino-sphere with different angles acquire different 
phases that decohere the neutrino ensemble, rendering less collective or even suppressing the neutrino modes. While the ’’single-angle” 
approach captures main features in many cases, it can sometimes modify or miss important aspects. 
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e.g., [I11I3II71II1I1I1I3IH33])- The coherent-state path formalism is used in Ref. [35] where the mean-field equations 
are shown to correspond to the stationary phase of the path integral for the many-body system. Ref. |36| derives 
mean-held equations and dispersion relations to describe neutrinos in media using relativistic Wigner functions. Ref. 
m has employed the BBGKY hierarchy which gives a rigorous theoretical approach to obtain mean-held equations 
and to reduce the full many-body description to an effective one-body description. In Ref. m an approach based on 
SU(2) algebras is used to put in relation the neutrino Hamiltonian with mixings and self-interaction with the (reduced) 
BCS Hamiltonian. This is extended to the three havor case in Ref. [38] . This approach shows that the corresponding 
full many-body probleirj^is exactly solvable. The many-body nature of the neutrino propagation problem in presence 
of self-interactions is also discussed in a schematic model based on SU(2) in Ref. [55] . 

Quantum kinetic equations are obtained in Refs. nanaiiaEi] based on a perturbative expansion of the interaction 
terms. Mean-held contributions of the usual typ^as well as the Bolzmann collision terms with the ’’molecular chaos” 
ansatz are included. A different derivation of these equations is presented in Ref. m where hrst-quantization and a 
perturbative expansion of the S matrix are used. In this case, the statistical factors of the Boltzmann term are not 
included and need to be put ” by hand”. The CTP formalism for out-of-equilibrium quantum helds has been used in 
Refs. [141 134] . Equations of motion for the neutrino Green’s functions are given, in the Boltzmann approximation. 
Ref. [10] calculates lowest order corrections to the dispersion relation and refractive index for neutrinos in a thermal 
bath. Tadpole and bubble diagrams contributions to the neutrino self-energy are calculated using lowest order thermal 
Green’s functions. In particular the effects are included of gauge boson masses and the imaginary part of the refractive 
index to account for neutrino scattering. It is shown that, in an almost GP symmetric bath, corrections from the gauge 
bosons masses dominate over the usual ones, while in supernovae the asymmetry between matter and anti-matter 
render them negligible compared to the usual mean-held corrections. 

On the other hand, current supernova simulations implement neutrino transport through the solution of the rela¬ 
tivistic Boltzmann equation 


m=c[f] ( 1 ) 

for the neutrino distribution functions / without taking into account mixings. In the most general case the distribution 
functions depend on seven variables t,r,p. The relativistic Liouville operator C is given by p^5f/5t + p’‘Sf /Sx'^ — 
T^y^yP^p''5f /5p^ with the Ghristoffel symbols, C the collision operator (see for example Refs. [6l I4T1I43] 1. Liouville 

equations for relativistic neutrino distribution matrices are derived in Ref. [20] . 

If one is interested in the mean-held approximation only, simple methods can be used to derive equations, as for 
example ttose in Refs. dZlElllSSI- We present the procedure and the main results obtained in Ref. m which has 
provides the most general mean-held equations. 

III. THE DENSITY MATRIX APPROACH 

In the mass basis, at each time the spatial Fourier decomposition of a Dirac neutrino held reads 

iP,{t,x)=[ e''P'^ij^{t,p,s), (2) 


with 


'ipt{t,p,s) = ai{t,p,s)ui{p,s) +bl{t,-p,s)vi{-p,s), 


(3) 


^ Note however that the neutrino full many-body Hamiltonian considered in the connection with BCS has mixings and neutrino self¬ 
interaction in the ’’single-angle” approximation. There is no matter term and the self-interaction coupling is constant|37|. 

® ’’Usual” means here without non-zero mass corrections and neutrino-antineutrino correlation terms. 
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where we note = J and Dirac spinors corresponding to mass eigenstates i are normalized 

as (no sum over i) 

uj (p, S)U^ {p, s') = vl {p, S)V, {p, s') = Sss' ■ (4) 

The standard particle and antiparticle annihilation operators (in the Heisenberg picture) for neutrinos of mass 
momentum p and helicity s satisfy the canonical equal-time anticommutation relation^ 

{ai{t,p,s),a]{t,p' ,s')} = {2Tr)^S’'^\p-p')Sss'Sij (5) 

{ai{t,p,s),aj{t,p',s')} = {al{t,p,s),a]{t,p',s')} = 0 (6) 

and similarly for the anti-particle operators. 

In the flavor basis, the field operator is obtained as 

X) = Uai 1pi{t, X), (7) 

with U the Maki-Nakagawa-Sakata-Pontecorvo unitary matrix [47]. Note that the indices can refer to active as well 
as to sterile neutrinos. In the framework of three active neutrinos the three mixing angles of U are now determined. 
Two are almost maximal, while the third one is small. The Dirac and Majorana CP violating phases are still unknown 

BHj. 


A. Two-point correlators 


Flavor evolution of a neutrino, or of an antineutrino, in a background can be determined using one-body density 
matrices, namely expectation values of bilinear products of creation and annihilation operatorj^ 

Pij{t,q,h,q',h') = {a]{t,q',h')ai{t,q,h)), ( 8 ) 

q, h, q', h') = {b] (t, q, h)bj{t, q', h')), (9) 

where the brackets denote quantum and statistical average over the medium through which neutrino propagate. For 


particles without mixings, only diagonal elements are necessary and relations (8]|9) correspond to the expectation 


values of the number operators. If particles have mixings as is the case for neutrinos, the off-diagonal contributions 
{i ^ j) of p and p account for the coherence among the mass eigenstates. 

The mean-field equations employed so far to investigate flavor evolution in astrophysical environments evolve the 
particle and anti-particle correlators p and p. However, the most general mean-field description includes further 
correlators. First, densities with ’’wrong” helicity states, such as pij{t,q,—,q',+) are present. These have already 
been shown to impact neutrino evolution in presence of magnetic fields [T5|. or if non-zero mass corrections are included 
[HISQ]. Moreover two-point correlators called abnormal or pairing densities [13 [HI, 


Ky (t, q, h, q', h') = {bj{t, q', h')a^{t, q, h)), 
Klnit,Q,h,q',h') = {aUt,q',h')bl{t,q,h)), 


( 10 ) 
( 11 ) 

also exist. If neutrinos are Majorana particles, correlators similar to ( To][TT ) can be defined, as done in Ref. such 
as {aj{t,q',—)ai{t,q,—)) or {b^j{t,q',+)bl{t,q,+)) that violate total lepton number. Equations of motion for these 
have been derived.in Ref. m Pair correlations are the fermionic analog of squeezed bosonic states (see also Ref. [H]). 


^ Note that a Fock space for flavor states can be built |44|. In the inflnite volume limit the flavor and mass operators belong to two 
representations that are not unitarily equivalent. In fact, although both satisfy the canonical anti-commutation relations, the vacua do 
not belong to the same Hilbert space [4^. This fact has been the object of debate (see, e.g., Ref. 061). 

A different convention pij{t, q, h, q', h') = (t, q, h)bi(t, q', h')) is also used in some works. With the one used here, the anti-neutrinos 

transform in the same way as neutrinos under the U transformation. 
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B. General mean-field equations 

The effective mean-field Hamiltonian takes the general bilinear form (h = c = 1) 


( 12 ) 


where ipt denotes the f-th component of the neutrino field in the mass basis Eq. ([^. The explicit expression of 
the kernel T depends on the kind of interaction considered (charged- or neutral-current interactions, non-standard 
interactions, effective coupling to magnetic fields, etc...). However it does not need to be specified to obtain the 
general structure of the equations. 

Equations of motion for the neutrino density matrix Eqs. (8]|9| can be obtained from the Ehrenfest theorem: 


ipij{t,q,h,q\h') = q', h')a^{t, q, h), H^sit)]) 

and similarly for the other correlators. Spinor products can be introduced 

r,7(t, q, h, q', h') = Ui{q, /i)fy (t, q- q')uj{q\h'), 

Q, h, q', h') = v,{q, h)fy (t, -q + q')vj{q', h'), 
rr/ (t, q, h, q', h') = Ui{q, h)f,j(t, q + q')vj{q', h'), 
rr/ (t, q, h, q', h') = v,{q, h)fy (t, -q- q')uj{q', h'), 


(13) 

(14) 


where the Eourier transform of the mean-held in Eqs. (14) is dehned as 

r,,(f,f)= f 

Jp 

It is straightforward to show that their evolution is determined through the following equations: 

zp{t) = r-'-'it) ■ Pit) - pit) ■ T-'-'it) + T-^'^it) ■ «t(t) - «(t). r‘^^(t), 

= v^^t) ■ Pit) - m ■ r""(i) - T'^'^it) ■ Kit) + K\t) ■ r‘^"(t), 

tkit) = r^'^it) ■ Kit) - Kit) ■ T^^it) - T'^'^it) ■ pit) - pit) ■ r'^^t) + r'^^it) 


(15) 


(16) 


An equation similar to the one for «;(t) holds for K^it). These are the most general mean-held equations for massive 
neutrinos prop aga ting in an inhomogeneous medium. 

Erom Eqs. (16) one can see that the evolution equation for k (and k^) depends on p and p as well as T'^'^ which 
implies that the usual consideration that the k correlators should oscillate fast around zero does not necessarily 
hold. Moreover, even if the k terms are zero initially, the particle-antiparticle mixing component of the mean-held 
Hamiltonian r'^^(t) acts as a source term producing non-zero k that can have a non-trivial back-reaction on p and p. 

Hollowing Refs. m and Eqs. (|16[) can be cast in a compact matrix form: 


tTZit) = [Hit),nit)], 


where the generalized Hamiltonian is 


and the generalized density 


nit) = 


nit) = 


Y^’^it) r^p(t) 
T^''it) T^^it) 


pit) Kit) 

K\t) I -pit) 


(17) 


(18) 


(19) 






Such equations become more explicit when one makes assumptions about the background. A common hypothesis, of 
interest for various applications, is the one of a homogeneous, (an)isotropic, unpolarised medium. Such a condition 
for the background corresponds to 


p{t,p'h',p,h) = {2Tr)^Shh'S^{p- f)pit,p), 


( 20 ) 


where p here corresponds to the particle composing the background like electrons or neutrinos (see Eq.(47l). Using 
Eqs. ( 12p9| one gets the following components for the generalised Hamiltonian H Eq.(18l for massless neutrinoj^ 


r’^‘'{t,q) = Sit,q)-q-Vit) 




( 22 ) 


with having the same expression as but with the V (t) contribution having a plus instead of a minus sign. The 
unit vectors q and Cq point to the neutrino direction of motion and perpendicular to it, respectively. The Nf x Nf 
scalaiH and vector matrices arj^ 

S{t, k) = h\q) + /i“"‘(t) + ^Gf f e{t,p). (23) 

Jp 


and 


V{t) = v^^\t) + V2GFj^\^pi{t,p) + epK{t,p) + e*K'l'(t,p)|, 


(24) 


where we defined 


p{t,-q). 

The scalar and vector matter contributions in the active neutrino sector read, in the flavor basis, 

1 


Ne(t)5^e - 


= V2GFdo,p 

U„7*(t) = V2GFd^i3 

with the particle number and velocity densities {vf = p/E^) of the particles composing the medium 

Nf{t) = 2 [pf{t,p) and J/(t) = 2 /i;/p/(<,p). 

J p J p 

The final expression for the mean-field Hamiltonian in its 2Nf x 2Nf matrix form thus reads 

f Sit,q)-q-Vit) -e*q-V{t) 

V S{t,q)+q-V{t) 


(25) 

(26) 

(27) 

(28) 

(29) 


One can see that the pairing correlations introduce neutrino-antineutrino mixing through the off-diagonal vector term 
which gives a contribution perpendicular to the neutrino momentum. This vector contribution is non-zero in presence 


The light-like four-vectors are defined as 

n^iP) = (^) ^ ( sp ) ’ 

where p = p/p denotes the unit vector in the direction of p and the pair of complex vectors (cp, e*) spans the plane orthogonal to p, with 
ip • €p = 0, €p • e* = 2. In terms of an oriented triad of real orthogonal unit vectors (p,for instance the standard unit vectors 
associated to p in spherical coordinates, one has ip = pg — ip^p. 

An extra scalar contribution is needed for the kernel in the Majorana case as pointed out in Ref. 09]. 

The expression for S is the same as for S but with a minus sign for the contribution. 
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of anisotropy. Note that contributions from neutrino-antineutrino pairing correlations have been first introduced in 
Ref. [in] using BBGKY (see Section 3.3). 

Under different approximations one recovers the mean-field equations usually employed in flavor studies in media. 
First, if neutrino-antineutrino correlations are neglected as done so far, the off-diagonal contributions to the total 
Hamiltonian % of Eq. (29) are zero, implying that there is no neutrino-antineutrino mixing. Note that an order-of- 


magnitude estimate for one neutrino generation indicates that the ratio of the neutrino-antineutrino mixing to the 
difference of the diagonal elements is likely to be small [49]. 

In the case of the Sun, the neutrino self-interaction contribution is negligible, the medium is homogeneous and 
isotropic with good approximation. One then finds from Eqs. (161, a Liouville Von-Neumann equations of motion for 
the neutrino and antineutrino density matrices: 


ip{t) = [h,p] ip{t) = [h,p]. 

with h = + h™®* and h = —hP + j^mat wherthe matter contribution is 

^mat ^ ^GpPe. 


(30) 


(31) 


which is the well known mean-held Hamiltonian that is responsible for the MSW effect |3]. In presence of the neutrino 
self-interaction Eqs. (30) receive a supplementary contribution 


Ki, = S{t, k) - q-V{t) = V2 Gf f {1 - q ■ p)i{t,p). 

J V 


(32) 


that has been shown to produce collective neutrino havor conversion modes iHHiniiiij- 

The inclusion of contributions from the neutrino mass is straightforward in the present approach (detailed derivations 
for Dirac and Majorana neutrinos are given in Ref. m- From Eqs. ( |16| ), by taking the homogeneity condition Eq. 
(20), the inclusion of the lowest order corrections from the neutrino mass (or corrections from the neutrino magnetic 
moment) requires to consider the helicity structure of the correlators (|8|^ as well 


p{tA) 


p—it,q) p-+it,q) \ ^ C{t,Q) 

P+-it,Q) P++it,Q), J ~ \ Kt,?) 


(33) 


where each quantity is to be understood as being a 2Nf x 2Nf matrix in mass/havor and helicity space. In the case 
of Majorana neutrinos, the total Hamiltonian including mass contributions is: 


where 




HM{t,q) = S{t,q) -q-V{t)-q- 
HM{t,q) = S{t,q) + q-V{t) + q- Vmit), 

m m 


<l>M(t,g) = e*^«e* 

with the supplementary off-diagonal contribution $m(^, ?) and off-diagonal contribution 




V^{t) = -V2GFj^{e-^^-epn{t,p)^+h.c.Y 


(34) 

(35) 

(36) 

(37) 

(38) 


Note that similar contributions are obtained for the neutral current on electrons, protons and neutrons. The contributions from electrons 
and protons exactly cancel because of the neutrality of the medium, while the ones from neutrons give rise to a term proportional to 
the unit matrix which does not affect flavor. 









10 


One can see from Eqs. ( Mp8 ) that in the ultrarelativistic limit, the evolution of the neutrino and anti-neutrino sectors 
decouple as expected. On the other hand, non-zero mass corrections introduce neutrino-antineutrino mixing. Since its 
contribution is also vectorial in nature, it requires an anisotropic medium to be non-zero. This is the spin coherence 
(also called helicity co/ierence[T7]) that has been pointed out in Ref. [T3] using the CTP formalism and the 2 Particle- 
Irreducible (2PI) effective action. A first schematic calculation with one neutrino flavor shows that corrections from 
the neutrino mass might produce flavor change through the presence of a MSW-like resonant conversion|5()j. On the 
other hand if one considers the neutrino magnetic moment (/r) coupling to stellar magnetic fields {B), the equations 
of motion acquire the same structure as in Eq. (34) but with the replacement = I^B. Ref. m has shown that 
significant neutrino-antineutrino flavor conversion can occur, even for a small magnetic moment compatible with the 
Standard Model, with nonzero masses and reasonable values of stellar magnetic fields. 


C. The application of BBGKY hierarchy 


A rigorous framework to go from the full many-body to the effective one-body description is given by BBGKY 
m- This is a hierarchy of coupled integro-differential equations for reduced density matrices or s-body correlation 
functions |52| that allows to derive from first principles the evolution equations at different levels of approximation, 
i.e. mean-field, extended mean-field, Boltzmann and beyond. These are obtained through different truncations of 
the hierarchy. Originally BBGKY was introduced for a N-body system of non-relativistic particles while it can be 
generalized to relativistic systems with an infinite number of degrees of freedom]^. In Ref. [TH], BBGKY has been 
applied for the first time to realistic neutrino systems. This furnishes a first principle derivation of the mean-field 
equations for neutrino matter and neutrino-neutrino interactions. Moreover the use of the BBGKY approach allows 
to go beyond. 

BBGKY is a hierarchy for s-reduced density matrices defined as {h = c = 1) 


= 


Nl 


{N-s)\ 


tr 


s+i...nD, 


(39) 


trs+i indicating that we are tracing over the s-|-I particle. D is the many-body density matrix that satisfies the 
Liouville Von-Neumann equation: 


i 


db 

dt 


[H,b]. 


(40) 


The Hamiltonian of the system H comprises both a free and a two-body interaction term, i.e., H = X)fc'^o(^) + 
J2k<k' b{k,k'). Instead of solving Eq. (1^), the exact evolution of the many-body system can be determined by 
solving the following hierarchy of integro-dinerential equations for the s-reduced densities: 

' ipi = [Hq{ 1), pi] + tr2[V{l,2), P 12 ] 
ipi2 = [Ho{l) + Ho{2) + V{l,2),pi2] 

+tr3[K(l,3) + I/(2,3),pi23] 

< . (41) 

*Pi...s = ELi^o(^) +Efc'>fe=i^(^>fc'),Pi...^] 

wherj^ e.g., 

Pi = (al“i)> P 12 = (4«i®i“2), (42) 


Similar definitions hold for the three-, four-, s-body reduced density matrices. 
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are the matrix elements components of the one-body and two-body reduced densities respectively (1, 2 indicate here 
particle 1, particle 2, etc...). Eqs. (41) represents an unclosed set of equations that couples the s- to the s + 1-system 
via the last term. It can be re-written as a hierarchy for connected correlation functions [52] or equal-time many-body 
Greens’ functions. For a system of relativistic particles the hierarchy becomes an infinite set of equations [S3]. 

Within BBGKY, the mean-field (Hartree or Hartree-Fock) approximation can be obtained by separating the cor¬ 
related and uncorrelated contributions to the two-body density matrix, namely, 


Pi2 — P1P2 + C12, 


(43) 


where the first term is the uncorrelated part while the correlated one is given by C 12 , the two-body correlation function. 
If the correlated part of the two-body density is neglected, then pi2 = P1P2 thus giving the following equation for the 
evolution of pi 


ipi = [Ho{l),pi] + tr 2 [V{1,2), PiP 2 ] = [hi{p),pi], (44) 

with the mean-field ri(p) = tr2{V{l,2)p2) produced by particles of type 2 acting on particle 1. More explicitly this 
is given by: 


ri.ij(p) = 


^)P2,■ 


(45) 


Expression ( |45| ) shows that the mean-held depends on the amplitude of the scattering process since the mean-held 
potential is built up from a complete set of one-body density matrix components for particle 2, p 2 ,nm, each contributing 
with the matrix elemen1p^U(im,jn) = (*w|Vi 2 |jn), with jn {im) incoming (outgoing) single-particle states. Going to 
the inhnite volume limit (uie discrete sum becomes an integral over the degrees of freedom of the particles constituting 
the background (expression (45) corresponds to ’’closing the loop”). As for the initial conditions, if one starts with 
an independent single particle state and neglect the correlated part of the two-body density Eq. (43), the many-body 
system keeps an independent particle system at all times. 


To go beyond the mean-held approximation contributions from the equation for pi 2 in (41), or correspondingly C 12 , 
need to be included. Such an equation reads: 


*Cl2 = [hl{p) +h2{p),Ci2] 

+ (1 - Pi){l - P2)V{1, 2)piP2{l - P 12 ) - (1 - Pl 2 )plP 2 V{l, 2)(1 - Pi){l - P 2 ) 
+ (1 -Pi- P2)V{1, 2)ci2 - Ci2E(1, 2)(1 -pi- P 2 ) 

+tr 3 [V{l, 3), (1 - Pi3)piC23(l - P 12 )] + trs[V{2, 3), (1 - ^ 23 )^ 2013(1 - P 12 )], 


(46) 


with Pij is the operator that exchange particle i with j. The hrst term on the r.h.s. corresponds to the mean-held 
acting on particles 1 and 2, the two terms on the second line account for collisions among particles, the third line 
gives contributions coming from correlations themselves and the last line corresponds to three-body terms with tr^ 
being the trace over the third particle. The common practice is to neglect the three-body terms since contributions 
from higher rank correlation functions are expected to be decreasingly importanip^ By retaining contributions from 
collisions only (the second line of (46)) and making the ’’molecular chaos” ansatz - the time between collisions is too 
short for the correlations between collisions to build-, one obtains a quantum Boltzmann equation for a system of 
neutrinos, consistently with, for example. Ref. m- 

Applying BBGKY to a realistic case for a system of neutrinos and anti-neutrinos requires the extension of Eqs. 
(41) to density matrices (or two-point correlators) p and p Eqs. ([8]|^. Using Eqs. ( 44][45 ), the mean-held associated 
to neutrino matter and neutrino-neutrino interactions is calculated. The low-energy limit of the Standard Model 


Note that in case of identical particles the matrix elements are antisymmetrised, i.e. = (im\Vi 2 \jn) — {im\Vi 2 \nj)- 

Note however that there are contexts in which the inclusion of three-body terms is necessary (see e.g. Ref. [53]). 
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FIG. 1. The figure shows, in a pictorial way, the interaction terms and the corresponding mean-fields associated with the 
neutrino-neutrino interaction. The two contributions correspond to the diagonal part of the mean-field arising from the usual 
scattering terms (upper figures) and the off-diagonal part (lower figures) associated with Pantaleone off-diagonal refractive 
index [19]. 


Lagrangian is sufficient for applications in the astrophysical context considered herq^ The neutrino self-interaction 
mean-field has both diagonal and off-diagonal contributions as first pointed out in Ref. [21] (Figure 1): 


= + k - p' ,h'J\ 

zV ^ Jp,p' 

{p, /ia) 7 '"(l - 75)^1.^ (t h'is)] (al^ (p, ha)a^f, {p', h!^)), 


(47) 


Mean-field equations for neutrinos have been derived by taking the homogeneity assumption (20) for the expectation 
value in Eq.(??). The quantum and statistical average over the full many-body system, in the definition of the 
s-reduced density matrices, reduce to expectation values over single particle states, giving the evolution equations 

(Isopl. 


If the medium in which neutrinos are propagating is dilute enough to neglect collisions, one can still have con¬ 
tributions to C 12 coming from third line of Eq. (461. Among po ssib le two-body cor relations are those of the type 
C 12 = K, ■ from neutrino-antineutrino pairing correlators k Eq. (10) and Eq. (11) that can give rise to a pairing 


mean-field. This can be defined similar to Eq. (45) with p being replaced by k, and the sum being over the initial 
(or final) single particle states [12]. The homogeneity condition for k is similaJ^ to the one for p Eq. (20). Such 
a condition implies that the neutrino and anti-neutrino operators correspond to particles with opposite momenta. 
Extended mean-field equati ons have been derived including such corrections for Dirac neutrinos whose structure is 
consistent with Eqs. ( IfiflO ), although derived for homogeneous and unpolarised media. Note however that in the V 
term Eq.(24) only the contribution from k and are non-zero while it has been pointed out in Ref. m that V is 
also sourced by the usual densities through the i terrrp^ 

Moreover it has been first shown in Ref. [55] that the generalized neutrino Hamiltonian with pairing correlations 
can be diagonalized by generalized Bogoliubov-Valatin transformation from the particle to the quasi-particle degrees 
of freedom. Such quasi-particle operators combine neutrino and anti-neutrino operators with opposite momenta. A 
special case for such transformations is represented by the special Bogoliubov-Valatin transformations: 


“fc = vial + ulbj, 
“I = ^kal + Wkh^ 


(48) 


Note that corrections from the gauge boson propagators need to be considered in the cosmological context |40|. 
Note that 6^{p — p') is replaced by 5^{p -\- p'). in this case. 

This difference is due to the fact that in Ref. m the number of particles has been kept fixed on average. 
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with k and k corresponding to opposite momenta. Similar relations hold for the hermitian conjugate operators. The 
anti-commutation rules imply that the Uk, Vk , Wk, Zk coefficients have to satisfy the following relations: + = 1, 

\w\^ + \z\^ = 1, |m|^ -I- \w\1 = 1, \z\l -I- \v\l = 1 and VkZk + WkUk = 0, ZkW^ + v^Uk = 0. In the limit of vanishing 
pairing mean field, since —)■ 1, Ufc —>■ 0 and Wk —>■ 1, Zfc —>■ 0, the quasi-particle operators and tend to the 
antineutrino operators and bj.; while the quasi-particle operators a\. and Uk tend to the neutrino operators a\ and 


a\. In this case, the extended Hamiltonian with pairing correlations reduces to the usual mean field. 

Determining neutrino propagation in astrophysical environments such as core-collapse supernovae is involved nu¬ 
merically, especially with the inclusion of neutrino self-interactions, of turbulence, of realistic geometries for the 
neutrino emission and of the background in which neutrinos propagate. The use of linearized methods provides a 
simple tool to investigate the occurrence of flavor collective modes and of flavor instabilities. A linearization of the 
neutrino equations of motion and of the eigenvalue equations are given in Ref. [56] . In particular one assumes that 
the system is in a quasi-stationary state and performs small amplitude variations around it in flavor space. It is in 
fact well known (see, e.g., Ref. m) that from Eqs.( |41[ ) linearised versions of the equations of motion can be obtained 
that are of great usa^ Ref. [5S] has given a different derivation of the linearized equations following methods known 
in the study of many-body systems such as atomic nuclei and metallic clusters. Such linearized eigenvalue equations 
can be cast in the well known form: 


[A B 
yB A 



(49) 


where the off-diagonal contributions to the neutrino p' and antineutrino p' density matrices correspond to the forward 
and backward amplitudes of the random-phase approximation. The presence of collective stable or unstable modes 
depend on the eigenvalues being real or complex. The matrix on the left hand side is the stability matrix that is 
related to the energy density curvature. Its concavity (or convexity) informs about the stability of the small amplitude 
collective modes. The A and B matrices depend in particular upon the derivatives of the mean-field Hamiltonians on 
the densities (see Ref. [55] for details). 


IV. A COHERENT-STATE PATH-INTEGRAL APPROACH 


An alternative way to obtain mean field equations exploits the path-integral formalism |35j . The mean-field equa¬ 
tions including matter and neutrino self-interactions for supernova neutrinos are shown to correspond to the saddle 
point approximation of the path-integral for the entire many-body system. To this aim the evolution operator can 
be determined by seeking a path integral representation for the neutrino Hamiltonian including the mixing terms as 
well as the neutrino interactions with matter {H^ in this section) and neutrino backgrounds An algebraic 

formulation of the neutrino Hamiltonian is performed in terms of the following generators of SU(2) algebras^ 


J+{p) = al{p)ae{p), J-{p) = al{p)ax{p), Jo{p) = ^ (4(p)aa;(p) - al{p)ae{p)) , 


(50) 


with aj.(p) {ax{p)) creation (annihilation) operator for a neutrino of flavor x = p,T and momentum p (similarly for 
the electron flavor e). The operators (50) satisfy the commutation relations 

[J+(j)),J-{q)]=2d^{p-q)Jo{p), [Jo{p), J±{q)] = ±S^{p - q)J±{p)- (51) 


In the context of many-body physics Eqs. | |30|32[ | are called Time-Dependent-Hartree-Fock (TDHF) equations. Their linearisation gives 
the Random-Phase-Approximation. When contributions from the pairing density k are included, one obtains the Time-Dependent- 
Hartree-Fock-Bogolioubov approximation. Their linearized eigenvalue equations are known as the quasi-particle random-phase approx¬ 
imation (QRPA). 

Here we restrict to the case of two flavors, while the symmetry of the three flavors case becomes manifest by using an SU(3) algebraic 
formulation m- 
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There are as many as the number of neutrino momenta in the situation of interest. Using the SU(2) coherent states 

\z{t)) = iVexp (52) 

with N a normalisation constant, V the ensemble of allowed momenta. The state \(f)) is 

!<(') = n (53) 

p&V 

|0) being the particle vacuum, one has the path-integral representation: 

{z\tf)\U\z{U)) = I D[z,z*]e^^^^’^‘\ (54) 

for the matrix element of the evolution operator between the initial (ti) and final (tf) times. The path-integral 
measur^EE] is 

D[z,z*] = hm ff TT 2 &iP^^o)dz*{p,t^) ^ 

^ N^oc J-J- J-J- 2m ^ ’ 

Q;=l 


The leading contribution to the path-integral 

S[z,z*] = - Hv- Huu) - i\og{z'{tf)\z{tf)) 


(56) 


comes from the stationary path |z(t)) that minimizes the action functional and can be determined solving the Euler- 
Lagrange equations. One obtains first-order non-linear equations of the Riccati-type for z{p,t). By interpreting it as 
ratio of one-body neutrino amplitudej^ 


z{p,t) 


'<PeiP,t) ’ 


(57) 


with the normalisation + IV’aiP = 1 expressing probability conservation, the equation for z{p,t) can be rewritten 
as the Schrodinger-like equation 


. d ( A + B — Is. cos 29 Bex + A sin 29 \ / tj;e\ 

^ dt y-ilJx ) ~ 2 \ B^e -I- A sin 20 -A-B + /Scos29j\i^:e)' 


where A = and A = \/2GFNe. B and Bex are given by 


B = 


y/2G[ 

V 


(fqRpq [|■^/'e((7,^)|^ - |V'x(g,0n ) 


(59) 


Bex — 


2V2Gf 

v~ 


cfqRpqil;e{q,t)i’x{(lA), 


(60) 


The exponential factor in the measure arises because the coherent states are over-complete. 

Here 'ipe{p,'t) and ipx{p-,t) are the electron neutrino disappearance and non-electron neutrino type appearance amplitudes, respectively. 
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with Rpq = (1 —cos'iJpq) the angular dependence coming from the spinorial products associated with the V — A nature 
of the weak interaction ('dpg being the angle between the two neutrino momenta), consistently with the angular part 
in Eq.(32). Eqs. ( 58p0| are the mean-field equations currently used for studies of neutrino flavor conversion in 
core-collapse supernovae, written in terms of neutrino amplitudes instead of density matrices as Eqs.( 30][32 ). (To 
see the equivalence among the theoretical treatments in terms of density matrices, neutrino amplitudes, isospin (or 
polarization) vectors, or evolution operators, see for example Ref. [bOj.l Eqs. ( 58p0 ) are consistent with previous 
results such as Ref. [58] . 

To go beyond the mean-field approximation, corrections to the saddle-point solution of the path-integral can be 
determined with a Taylor expansion of the action around the stationary path. This procedure leads to corrections 
to the matrix elements of the evolution operator in the form of formal determinants (Eqs. (58-60) of Ref. [35]) that 
depends on the second derivative of the action along the classical path. Such corrections correspond to the small 
amplitude (or RPA) formulation in this formalism. Note that the extended mean-field equations provided in Ref. |35j 
do not include corrections from the neutrino mass, neither from neutrino pairing correlations. The latter is due to 
the chosen specific form of the coherent state Eq.(52). 

The algebraic formulation of the problem of neutrino propagation in media obtained in Ref. |35| is further developed 
in Ref. [37] where the neutrino Hamiltonian with self-interaction (in the ’’single-angle” approximation, without the 
matter term, neutrino density fixed) is shown to have the same form as the reduced Bardeen-Cooper-Schrieffer (BCS) 
Hamiltonian for superconductivity |59| . The algebraic approach reveals the invariants and the exact solvability of the 
corresponding collective many-body and one-body neutrino Hamiltonian. The spectral split phenomenon is interpreted 
as a transition from a quasi-particle to a particle description. Based on SU(3) the approach is extended to three flavors 
in Ref. |38| and the exact solvability of the full many-body problem (within the same approximations as before) is 
demonstrated. Moreover the question of CP violation effects in supernovae is addressed. The existence of such effects 
is established in Ref. [51] while they are found quantitatively to be small |STJ[B2|. The condition for their presence can 
be understood in terms of a factorisation condition of the neutrino Hamiltonian, as shown in Ref. m and extended 
in Ref. [62] in presence of neutrino self-interactions. CP violating effects can arise when the factorisation condition is 
broken, for example in presence of radiative corrections, of non-standard interactions, or of magnetic moments. Ref. 
[38] has shown that such a factorisation condition holds also for the many-body Hamiltonian (and not only in the 
effective mean-field description). 


V. THE QUANTUM BOLTZMANN EQUATION FOR NEUTRINOS 

An early derivation of the Boltzmann equation for relativistic distribution functions with mixings is made in Ref. 
m- Equations of motion in terms of density matrices for neutrinos and antineutrinos propagating in a background 
are first given in Ref. m- Collisions of a(n) (anti)neutrino with the particles composing the background are taken 
into account in a perturbative approach. These results are extended in Ref. m where the explicit expression for 
the diagonal and off-diagonal contributions of the collision integral are given both for charged and neutral current 
neutrino interactions with the medium. The derivation relies on the ’’molecular chaos” ansatz. This gives Boltzmann 
equations for particle with mixings 1 1 3] : 

ip{t) = [h,p] -f C[p,p] ip{t) = [h,p] -f C[p,p], (61) 

where the (7, C terms introduce collisions that have the role of bringing the system to thermodynamic and flavor 
equilibrium. These include in particular scattering on neutron and protons, electrons and positrons as well as pair 
annihilations since the equations are derived for neutrino evolving in the early universe plasma. Sometimes a more 
economical damping approximation is used where the ansatz is made that collisions drive the system to equilibrium 
exponentially (see for example Ref. m)- 

The quantum Boltzmann equations obtained in Ref. |13j have also been formulated using (iso)spins related to the 
flavor amplitudes. They give a useful way to picture the evolution in flavor space. To this aim, the Schrodinger-like 
equation for the neutrino evolution is replaced by a precession equation for the flavor isospins subject to effective 
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magnetic fields [64) . In the context of neutrino physics the approach is first applied in the early universe context 
m, then to the MSW effect [SS] and recently to picture flavor conversion phenomena in presence of the neutrino 
self-interactions (see Refs. [8l[24] and references therein). 

A general formalism to rigorously derive flavored quantum Boltzmann equations for isotropic systems is presented in 
Refs. based on the closed-time path (CTP) or ”in-in” formalism and the 2PI effective action (see Ref. [50] for 

a review). The approach has been applied to a model system of fermionic and scalar fields with a Yukawa interaction 
in the context of baryogenesis and leptogenesis. Contributions from neutrino-antineutrino pairing correlations are 
propertly taken into account and are shown to produce sizeable effects |28j . 

The application of CTP and the 2PI effective action to supernova neutrinos is performed in Refs. [nisi]. Equations 

of motion for the two-point functions, such as G'^‘}j{x,y) = (Tp(ipf{x)'tjjj°'{y))) are obtained for Majorana 


tao; 

neutrinos. Tp is the time ordering operator along a specific path. The CTP formalism considers a closed path 
starting at initial time to ap to the time of interest, and then back to the initial time. The 2PI effective action 
contains a one-loop contribution and a second contribution coming from all higher-loop terms. Only connected 2PI 
diagrams are included in the computation of the self-energy E in the Schwinger-Dyson equation: 


{i - M)G (x, y) - * y G y) = 1 [x - y), 


(62) 


where M are spin x flavor matrices. The two-point function can be decomposed into a spectral function, that encodes 
information on the particle states, and a statistical function, which carries information on the occupation numbers 
of these states. The spectral function varies little from the massless free held, while the dynamics of the statistical 
function is the one that interests us. In particular, two contributions to the Wigner transform of the statistical 
function are retainecj^ i.e., those from expectation values for the number operators for neutrinos and antineutrinos 


Its contribution is non-zero in presence of 


and from (^d}j (^i) bj {q 2 )J, the neutrino-antineutrino correlation functior 

anisotropy and gives rise to the spin coherence between the neutrino andantineutrino sector (see Section 3.2). While 
the Wigner transform of the equations of motion for the statistical function can be written as a gradient expansion with 
inhnite series of derivatives, these are truncated at lowest nontrivial order in gradients. The parameter e governing 
this expansion is actually the ratio of the neutrino wavelength over the typical (space or time) inhomogeneity. As 
a consequence, in the self-energies only one-loop contributions of the order of e appear in the local piece of the self¬ 
energy, while two-loop contributions (of the order of 0{e^)) contribute to the spectral and statistical components of 
the self-energy. Within these approximations equations of motion are derived for the neutrino densities matrices and 
spin coherence densities. The corresponding kinetic equations can again be cast in a compact matrix form: 

iD[T]-[n,T]=iC[J^]. (63) 


Here, for 3 neutrino flavors, J- and H are 6x6 matrices having the following block structure: 


F-( ^ ^ 

I fT 


n = 


H 

hU -h^ 


G = 


G G4, 
Gi 


The H^i, is given by 


= ~T (S+m* -I- . 


(64) 


(65) 


where one can see the pre sence of the vv mixing terms, consistently with Eq. (29) of Ref. [m, suppressed by m/k as 
expected (see also Eq.(35) for ^m)- The explicit expression for H and the derivative term D and the collision 
terms C, G and G^ can be found in Ref. mi. Only the vv —>■ w contributions to the collision term are shown. Given 
the approximation made in the calculation of the collision integral, the quantum Boltzmann equations obtained in 
Ref. mi with the full collision term should be consistent with the ones in Ref. m| based on the density matrices. 


Contributions from pairing correlations are neglected because they are expected to vary on small timescales. 

This kind of contribution is first considered in Ref. m in relation with the neutrino magnetic moment in presence of stellar magnetic 
fields. 
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VI. CONCLUSIONS 

Understanding how neutrinos change their flavor under different astrophysical conditions remains a fascinating 
problem. Flavor conversion in the Sun and in the Earth is well understood in terms of the MSW effect (also of 
the parametric resonance[571 [55] for the Earth). The case of core-collapse supernovae and of binary systems keeps 
revealing novel mechanisms. A variety of theoretical approaches have been used in the past two decades to describe 
neutrino propagation in such dense environments as well as in the early universe. 

In this review we have summarised the methods employed, with a particular emphasis on the commonly used mean- 
field equations including the neutrino mixing, neutrino-matter, neutrino-neutrino interactions and on the Boltzmann 
equations that include incoherent scattering. We have in particular highlighted features related to the many-body 
aspect of the problem as well as approaches in which such evolution equations can be derived from first principles. 
Obviously the methods can be employed to treat more general cases involving beyond the Standard Model physics such 
as non-standard interactions or novel particles. Several interesting developments are ongoing that have unravelled 
the connection between flavor conversion of neutrinos in an astrophysical background and properties of many-body 
systems such as atomic nuclei or in condensed matter. 

Since the equations of motion are non-linear, particular care is needed to verify the role of corrections to the mean- 
field. The study of the latter is necessary to put our description on solid grounds. In particular, pairing correlations 
(among like-particles or not) and spin (or helicity) coherence naturally appear in the most general mean-field de¬ 
scription of massive neutrinos propagating an anisotropic background. Further numerical calculations implementing 
helicity coherence should assess their actual role in realistic astrophysical conditions. The issue of pairing correla¬ 
tions might require relaxing the homogeneity condition or the inclusion of collisions to meet appropriate conditions. 
We emphasize that the equations for the pairing correlators are sourced by the usual neutrino density matrices and 
not only by the pairing correlations themselves. The inclusion of collisions in a realistic Boltzmann treatment with 
mixings is also a necessary step to investigate the competition of the three timescales of the problem, i.e., the time 
between collisions that introduce decoherence of the flavor evolution, the typical timescales for flavor change and the 
timescales over which the macroscopic system evolves. 

Moreover, an improved description, also of the transition region from the high to the low density regime in the 
supernova, requires a more realistic description of the neutrino emission at the neutrino-sphere, realistic geometries 
and inclusion of inhomogeneities in realistic conditions. The role of multi-dimensions, of convection and turbulence 
is clearly established in supernovae simulations. Linearized studies and schematic models of flavor evolution bring 
insight, for example in the search for instabilities in the heating region behind the shock, or of the possible role of flavor 
evolution on the supernova dynamics. However, the path ahead might definitely require realistic numerical simulations 
and extended theoretical descriptions. Such studies as well as the theoretical connections between neutrino flavour 
evolution in media, and other domains, might unravel interesting new aspects. 
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